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PREFACt 


This report describes EQLBRM, an Interactive F0R1RAN IV computer program 
for computing adiabatic flame temperatures and equilibrium product composition 
for the combustion of hydrocarbon fuels In air. The authors are grateful to 
Paul S. Cramer of the AVCO-Lycomlng Corporation for extensively testing ElJLBRM 
and suggesting many Improvements. 
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1. INTRODUCTION 


T .1 Overview 

EQLBRM Is a user-friendly, menu-driven. Interactive computer program which 
calculates the adiabatic equilibrium temperature and product composition 
resulting from the combustion of hydrocarbon fuels with air, at specified 
constant pressure and enthalpy. EQLBRM was developed primarily as an Instruc 
tlonal tool to be run on small computers to allow the user to economically and 
efficiently explore the effects of varying fuel type, air fuel ratio. Inlet 
air and/or fuel temperature, and operating pressure on the performance of 
continuous combustion devices Su, 1 * as gas turbine combustors, Stirling engine 
burners, and power-generation furnaces. 

1 .2 Method of Solution 

EQLBRM Is based on the principle of minimization of the Gibbs function of 
a reactive mixture at constant pressure, subject to constraints of conservation 
of atomic species and enthalpy. The minimization algorithm used Is the very 
efficient and robust Gordon and McBride method (ref. 1), as refined and sim- 
plified by Pratt (refs. 2 to 4), and recently developed Independently by 
Reynolds (ref. 5). A detailed description of the solution algorithm Is con- 
tained In reference 4. 

1 .3 Special Features 

The many "user-friendly" .'.'atives offered by EQLBRM Include: 

1. A resident default combustion problem specification allows Immediate 
demonstrc Mon of program execution and output format. 

2. Interactive, run-time selection of SI or English units. 

3. Fuel compounds may be either chosen from a resident data file, 
EQLBRM.DAT, or entered Interactively at run time. 

4. Composition of the adiabatic equilibrium combustion products Is out 
put In mole numbers, mole fractions, mass fractions and emission Index. 

5. At the user's option, a summary of the Inlet conditions together with 
the equilibrium solution may be written to a file, EQLOUT.LST, for subsequent 
hard copy output. 

1.4 Limitations 

Compared to the powe f ul NASA equilibrium code CEC71 ( 1 ) , EQLBRM has the 
following limitations: 

1. To fit In a small computer, program size Is limited to 80 species, 
compared to over 500 species In CEC71. (The limitation to 80 species Includes 
any new species added during an Interactive session.) 

2. Condensed species are not considered. 
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3. Only the HP problem (assigned enthalpy and pressure) Is considered by 
EQLBRM, while CEC71 considers also the TP, UV, TV, SP and SV problems (ref. 1). 

4 . Multiple problem statements are not considered. 

5. EQLBRM considers only the first NPROO species In the THERMO data file 
(see Chapter 3) as combustion products; NPROO Is presently hard-coded as 11. 
NPROD can be modified only by altering a DATA statement In the BLOCK DATA sub- 
program source code. 


2. USER'S GUIDE 

A well-written Interactive program contains the equivalent of a User's 
Guide In Its Interactive menu structure and prompts. As a practical matter. 

It Is necessary to actually use on Interactive program to discover how It 
works. One way to resolve this "Catch-22'' dilemma Is to obtain telephone- 1 Ine 
access to a demonstration run on a computer In which the Interactive program 
Is Installed. By this method a prospective user can "take a test drive" and 
see whether or not the program Is of further interest. 

An alternative method Is the presentation of example Interactive sessions. 
Two s'ch example sessions are presented heie. In both examples, everything 
that would appear on the user's terminal (either a printing terminal or 
cathode-ray tube screen) appears In UPPER CASE text, as It does In actual pro- 
gram execution. The user response to prompts Is boxed, and the lower case 
text are comments provided to rationalize the hypothetical user's choices, or 
to highlight some of the features of EQLBRM. 

2.1 Example Run 1 

In this example, It Is assumed that the user Is running EQLBRM for the 
first time, and wishes only to become familiar with Its operation, rather than 
solve a specific problem. The user therefore types In "RUN EQLBRM" or an 
equivalent system command to execute the compiled and linked load module named 
"EQLBRM." 

Immediately EQLBRM Introduces Itself by printing to the terminal: 

*** PROGRAM EQLBRM *** 

CALCULATES THE ADIABATIC EQUILIBRIUM TEMPERATURE AND PRODUCT COMPOSITION 

RESULTING FROM COMBUSTION OF HYDROCARBON FUELS AND AIR AT CONSTANT PRES- 
SURE AND ENTHALPY. 

USER MUST SPECIFY INITIAL FUEL COMPOSITION AND TEMPERATURE, INITIAL AIR 

TEMPERATURE, AIR-FUEL MIXTURE RATIO, AND PRESSURE. A VARIETY OF UNITS ARE 

AVAILABLE: 


TEMPERATURE 

KELVIN 

CELSIUS 

FAHRENHEIT 

RANKINE 


PRESSURE 
KILO PASCALS 
P.S.I. 
ATMOSPHERES 
BAR 


AIR-FUEL RATIO 
KG AIR/KG FUEL 
KG FUEL/KG AIR 
MOLS AIR/MOLS FUEL 
MOLS FUEL/MOLS AIR 
FUEL-AIR EQUIVALENCE RATIO 


v 
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EQLBRM WILL RETURN THE EQUILIBRIUM TEMPERATURE AND FUEL-AIR EQUIVALENCE 
RATIO, AS WELL AS THE MOLE NUMBER, MOLE FRACTION, MASS F RAC I ION AND EMIS 
SION INDEX OF THE PRODUCTS AT EQUILIBRIUM. 

PRESS <RETURN> TC CONT INUE 

Being an Interactive program, EQLBRM always prompts the user for the proper 
response or range of responses. The user responds »1th: 


freturrfr 


IF THE USER SPECIFIES A FUEL TOR WHICH EQLBRM HAS NO 1 HERMOCHE^ICAL DATA, 
THE USER MAY EXPLICITLY SUPPLY THE FUEL CHEMICAL FORMULA AND STANDARD 
STATE ENTHALPY OF FORMATION OR LOWER HEATING VALUE. AGAIN, A CHOICE OF 
UNITS WILL BE AVAILABLE TO THE USER. 

WOULD YOU LIKE TO SEE A SAMPLE DATA SET (Y/N) 

Since this Is just a familiarization run, the user elects to Inspect the 
sample data set. 


m 

‘“CURRENT DATA STATUS* * 


FUELS AND FUEL AMOUNTS : 

FUEL RELATIVE MOLES AMONG FUELS 

CO 1.0000000 

H2 2.0000000 

FUEL TEMP = 2.9800E+02 DEG- KELVIN 

AIR TEMP - 7 . 5000E+02 DEG- KELVIN 

PRESSURE = 5 . 0000E+00 ATMOSPHERES 

A-F RATIO = 8 . 5000E- 01 FUEL- AIR EQUIVALENCE RATIO 


SELECT AN OPTION: 

0- RUN WITH THIS DATA 

1- CHANGE FUELS OR FUEL AMOUNTS 

2- CHANGE FUEL TEMPERATURE 

3- CHANGE AIR TE"PEKATURE 

4- CHANGE TEMPERATURE UNITS 

b- CHANGE PRESSURE OR PRESSURE UNITS 
6- CHANGE AIR-FUEL RATIO OR RATIO UNITS 


OPTION NUMBER? (0-6) 


Again, since this Is Just a familiarization run, the user chooses to run with 
the sample data set. Just to see what the output looks like. 

0 

RUN WITM this DATA: ARF. YOU SURE? (Y/N) 

EQLBRM always double-checks before execution to avoid possible tine-costly 
user error. 
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Now EQLBRM executes the sample data set and returns the solution to the 
screen: 


****** 

OUTPUT 


*** ADIABATIC EQUILIBRIUM AT CONSTANT PRESSURE AND ENTHALPY *** 

PRESSURE= 5.000E+00 ATMOSPHERES 
TEMPERATURE = 2444.76 DEGREES KELVIN 

FUEL-AIR EQUIVALENCE RATIO = 0.850 


SPECIES 

MOLE 

MOLE 

MASS 

EMISSION 

NAME 

NUMBER 

FRACTION 

FRACTION 

INDEX 

("I") 

( MOL- I/G-X) 

(MOL-I/MOL-X) 

(G-I/G-X) 

( G-I/KG-FUEL) 

CO 

2.204E-04 

6.078E-03 

6 . 1 73E-03 

6 . 1 7 3E ♦CO 

C02 

3.270E-03 

9.019E-02 

1 . 439E-01 

1 .439E+02 

H 

1 .458E-05 

4.020E-04 

1 .469E-05 

1 .469E-02 

H2 

7.696E-05 

2.123E-03 

1 .551E-04 

1 .551E-01 

H20 

6.771E-03 

1 .868E-01 

1 .220E-01 

1 .220E*02 

N 

2.282E-09 

6 . 292E-08 

3 . 1 96E-08 

3 . 1 96E-05 

NO 

2.848E-04 

7.854E-03 

8.545E-03 

8 . 54 5E *-00 

N2 

2.417E-02 

6.667E-01 

6.772E-01 

6.772E+02 

0 

3.161E-05 

8.719E-04 

5 . 058E-04 

5.058E-01 

OH 

2. 503E-04 

6.903E-03 

4 . 257E-03 

4 . 257E+00 

02 

1 .163E-03 

3.209E-02 

3.723E-02 

3.723E*01 


DO YOU WANT A HARD COPY OF YOUR INPUT ANO OUTPUT? (Y/N) 

The user elects to make a hardcopy of both the Input and output datasets. 

0 

DATA AND RESULTS ARE BEING APPENDED TO 
FILE "EQLOUT . LST M IN YOUR DIRECTORY. 

The hard copy output Is directed to logical unit 2. System- specif 1c commands 
must be executed at run time (or added to the source code) to assign this unit 
to the file "EQLOUT. LST, " or directly to a printer. 

DO YOU WANT ANOTHER RUN? (Y/N) 


A "Y" response would return the user to the main Input menu, where the current 
Input conditions can be modified Interactively, as demonstrated In the next 
example run. However, the user in this case chooses to end the familiar- 
ization session 




2.2 Example Run 2 

In this example, we assume the user wishes to solve the following prob 
lem: Propane at 125° C burn* with air at 375° C, at constant pressure of 4.5 
atmospheres, with a 'lightly rich fuel-air equivalence ratio of 1.1. 

The user Initiates EQLBRM, and the program Introduces Itself: 

*** PROGRAM EQLBRM *** 

CALCULATES THE ADIABATIC EQUILIBRIUM TEMPERATURE AND PRODUCT COMPOSITION 
RESULTING FROM COMBUSTION OF HYDROCARBON FUELS AND AIR AT CONSTANT PRES- 
SURE AND ENTHALPY. 


USER MUST SPECIFY INITIAL FUEL COMPOS. TION AND TEMPERATURE, INITIAL AIR 
TEMPERATURE, AIR-FUEL MIXTURE RATIO, AND PRESSURE. A VAR1E1Y OF UNITS ARE 
AVAILABLE: 


TEMPERATURE 

KELVIN 

CELSIUS 

FAHRENHEIT 

RANKINE 


PRESSURE 
KILO PASCALS 
P.S.I. 
ATMOSPHERES 
BAR 


AIR-FUEL RATIO 

KG AIR/KG FUEL 

KG FUFL/KG AIR 

MOLS AIR/MOLS FUEL 

MOLS FUEL/MOLS AIR 

FUEL- AIR EQUIVALENCE RATIO 


EQLBRM WILL RETURN THE EQUILIBRIUM TEMPERATURE AND FUEL- AIR EQUIVALENCE 
RATIO, AS WELL AS THE MOLE NUMBER, MOLE FRACTION. MASS FRACTION AND EMIS- 
SION INDEX OF THE PRODUCTS AT EQUILIBRIUM. 


PRESS <RETURN> TO CONTINUE. 


(return) 


IF THE USER SPECIFIES A FUEL FOR WHICH EQLBRM HAS NO THERMOCHEMICAL DATA, 
THE USER MAY EXPLICITLY SUPPLY THF. FUEL CHEMICAL FORMULA AND STANDARD 
STATE ENTHALPY OF FORMATION OR LOWER HEATING VALUE. AGAIN, A CHOICE OF 
UNITS WILL BE AVAILABLE TO THE USER. 


WOULD YOU LIKE TO SEE A SAMPLE DATA SET? (Y/N) 


Since the user has a data set In mind already, there Is no Interest In seeing 
the sample dataset. 


0 

***CURRENT DATA STATUS*** 


FUELS AND FUEL AMOUNTS : 


FUEL 


RELATIVE MOLES AMONG FUELS 


- NO FUEL SPECIFIED 
FUEL TEMP * O.OOOOE-Ol 


AIR TEMP * 
PRESSURE = 
A-F RATIO = 


O.OOOOE-Ol 

O.OOOOE-Ol 

O.OOOOE-Ol 


DEG- KELVIN 
DEG- KELVIN 
ATMOSPHERES 

FUEL-AIR EQUIVALENCE RA1 10 
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SELECT AN OPTION: 

0- RUN WITH THIS DATA 
T- CHANGE FUELS OR FUEL AMOUNTS 

2- CHANGE FUEL TEMPERATURE 

3- CHANGE All TEMPERATURE 

4- CHANGE TEMPERATURE UNITS 

5- CHANGE PRESSURE OR PRESSURE UNITS 

6- CHANGE AIR-FUEL RATIO OR RATIO UNITS 

OPTION NUMBER? (0-6) 

Note that the main Input menu comes up blank. Data must be entered by the user 
Interactively. The order In which this Is done Is not Important, only It Is 
Important that no null data be left when the dataset is executed. For the 
sake of convenience, the user decides to begin with option 1, and work down: 

□ 

•CURRENT FUELS DATA* 

FUEL RELATIVE MOLES AMONG FUELS 

- NO FUEL SPECIFIED - 

SELECT AN OPTION: 

0- NO MORE CHANGES HERE 

1- CHANGE RELATIVE MOLES OF AN EXISTING FUEL 

2- ADD A NEW FUEL 

3- REMOVE AN EXISTING FUEL 

4- LIST ALL COMPOUNDS AVAILABLE IN THE THERMO LIBRARY 

OPTION? (0-4) 

The user Is now In the fuels sub-menu, and exercises the 2nd option - 
logically enough. 


0 

ENTER NEW FUEL NAME. 

The user could enter any alphanumeric string here - but EQLBRM Is looking for 
a chemical formula. We assume the user carelessly enters "C7H8" as propane: 

| C 7 H8 1 

ENTER RELATIVE MOLES OF C7H8 

If there Is to be a mixture of fuels, It Is necessary to specify the relative 
amount of each fuel. For the case of a single fuel, however, any nonzero 
value will do. The user chooses "1.0." 

1 . 00 | 

•CURRENT FUELS DATA* 
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FUEL RELAlIVt MOLES AMONG FUELS 

C7H8 T.OOOOE+OO 

SELECT AN OPTION: 

0- NO MORE CHANGES HERE 

T- CHANGE RELATIVE MOLES OF AN EXISTING FUEL 

2- ADD A NEW FUEL 

3- REMOVE AN EXISTING FUEL 

4- LIST ALL COMPOUNDS AVAILABLE IN THE THERMO LIBRARY 

OPTION. (0-4) 

The user notices the error that propane Is C3H8, not C7H8 as entered. Option 
3 must be exercised to remove the erroneous fuel se.ection. 

m 

DELETE C/H8 

-- ARE YOU SURE? (Y/N) 

EQLBRM double checks *1th the user In order to avoid time-consuming deletion 
errors . 


0 

‘CURRENT FUELS DATA* 

FUEL RELATIVE MOLES AMONG FUELS 

- NO FUEL SPECIFIED - 

SELECT AN OPTION: 

0- NO MORE CHANGES HERE 

1- CHANGE RELATIVE MOLES OF AN EXISTING FUEL 

2- ADD A NEW FUEL 

3- REMOVE AN EXISTING FUEL 

4- LIST ALL COMPOUNDS AVAILABLE IN THE THERMO LIBRARY 

OPTION? (0-4) 


Option 2 Is used again to enter the correct tuel name. 

0 

ENTER NEW FUEL NAMt . 


C3H8 


-- SORRY -- C3H8 

IS NOT IN THE THERMOCHEMICAL LIBRARY. IF YOU WISH TO USE 1 HIS COMPOUND 
YOU WILL HAVE TO SPECIFY THE ENTHALPY OF FORMATION OR LOWER HEAT L.LUE . 
DO YOU HAVE THIS INFORMATION (Y/N) 
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Thermochemical data for C3HE was not Included In che Input data deck. We here 
assume that the user has In mind a reasonable figure for the LHV of typical 
hydrocarbon. 


0 

EQLBRM Interrogates the user for the molecular composition of C3H8: 

SPECIFY CHEMICAL COMrOSlTION: 

ENTER KG-ATOMS C PER KG-MOL C3H8 


3.00 


EN'ER KG-ATOMS H PER KG-MOL C3H8 


8.00 


ENTER KG ATOMS 0 PER KG-MOL C3H8 

There Is no oxygen In the C3H8 molecule, so a null response Is required. 

lo.ool 


ENTER KG-ATOMS N PER KG-MOL C3H8 


0.00 


YOU MAY SPECIFY EITHER: 

1- STANDARD STATE ENTHALPY OF FORMATION 

2- LOWER HEATING VALUE 

WhICH VALUE? (1/2) 

The user recalls that a typ.cal LHV for hydrocarbon fueis Is about 18,500 
Btu/lbm, and decides to use this value. 


0 

SELECT UNITS: 

1- KJ/KMOL 

2- KJ/KG 

3- BTU/LB 

4- CAL/G 

WHICH UNITS? (1-4) 

0 

ENTER LOWER HEATiNG MLUE 


18500. 


ENTER RELATIVE MOLES OF C3H8 
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Again, since this Is the only fuel to be Included, any no.izero value will do. 


tun 

•CURRt.Nl FUELS DATA* 

F JEL RELATIVE. MOLES AMONG FUELS 

C3H8 1.0000E»00 

SELECT AN OPTION: 

0- NO MORE CHANGES HERE 

1- CHANGE RELATIVE MOLES OF AN EXISTING FUEL 

2- ADD A NEW FUEL 

3- REMOVE AN EXIS1ING FUEL 

4- LIST ALL COMPOUNDS AVAILABLE IN THE THERMO LIBRARY 


OPTION? (0-4) 

No more changes are required, but the user decides tc see what compounds are 
In the properties library: 

0 

•COMPOUNDS AVAILABLE IN THERMO LIBRARY* 


CO 

C02 

H 

H2 

H20 

N 

NO 

N2 

0 

OH 

02 

N02 

CH4 

C6H6 

C6H6( L) 

C7H8 

C8H16 

C8H18 

C8H1 8( L ) 

C12H26 

CH40 

C2H50H 

C3H80 

CH2 

CII20 

C3H8 

C2H2 

C2H4 

C2H6 


PRESS <RETU?N> TO CONTINUE. 

Note that C3H8 has been added to the list of recoqMsed compounds. It will, 
however, be lost when EQLBRM is exited. The user must add thermochemical data 
for C3H8 to the Input deck If C3H8 Is to be used repeatedly. EQLBRM can ac 
commodate up to 80 species. Inclusive of those added Interactively In the 
manner Just demonstrated. 


(return) 


•CURRENT FUELS OATA* 


FUEL RELATIVE MOLES AMONG FUELS 

C3H8 1 .OOOOEtOO 

SELECT AN OPTION: 

0- NO MORE CHANGES lit RE 

1- CHANGE RELATIVE MOLES OF AN EXISTING FUEL 

2 - ADD A NEW FUEL 

3- REMOVE AN EXISTING ‘UEL 

4- LIST ALL COMPOUNDS AVAILABLE IN THE THERMO LIBRARY 



OPTION? (0-4) 


No more changes are required. 


GO 

♦••CURRENT OATA STATUS*** 


FUELS AND 
FUEL 
C3H8 

FUEL TEMP 
AIR TEMP 
PRESSURE 
A-F RATIO 


FUEL AMOUNTS : 

RELATIVE MOLES AMONG FUELS 
1 .OOOOOOC 

- O.OOOOE-Ol DEG- KELVIN 

- O.OOOOE-Ol OEG- KELVIN 

- O.OOOOE-Ol ATMOSPHERES 

• O.OOOOE-Ol FUEL-AIR EQUIVALENCE RATIO 


SELECT AN OPTION: 

0- RUN WITH THIS OATA 

1- CHANGE FUELS OR FUEL AMOUNTS 

2- CHANGE FUEL TEMPERATURE 

3- CHANGE AIR TEMPERATURE 

4- CHANGE TEMPERATURE UNITS 

5- CHANGE PRESSURE OR PRESSURE UNITS 

6- CHANGE AIR-FUEL RATIO OR RATIO UNITS 


OPTION NUMBER? (0-6) 


The user wishes to work with Celsius as temperature units. 

□ 

SELELT AN OPTION: 

1- CHANGE UNITS TO KELVIN 

2- CHANGE UNITS TO CELSIUS 

3- CHANGE UNITS TO FAHRENHEIT 

4- CHANGE UNITS TO RANKINE 

OPTION? (1-4) 

m 

•••CURRENT OATA STATUS*** 


FUELS ANO 
FUEL 
C3H8 

FUEL TEMP 
AIR TEMP 
PRESSURE 
A-F RATIO 


FUEL AMOUNTS : 

RELATIVE MOLES AMONG FUELS 
1 .0000000 

- O.OOOOE-Ol DEG- CELSIUS 

- O.OOOOE-Ol DEG- CELSIUS 
* O.OOOOE-Ol ATMOSPHERES 

- O.OOOOE-Ol FUEL-AIR EQUIVALENCE RATIO 


SELECT AN OPTION: 

0- RUN WITH THIS DATA 

1- CHANGE FUELS OR FUEL AMOUNTS 

2- CHANGE FUEL TEMPERATURE 

3- CHANGE AIR TEMPERATURE 

4- CHANGE TEMPERATURE UNITS 

5- CHANGE PRESSURE OR PRESSURE UNITS 

6- CHANGE AIR-FUEL RATIO OR RATIO UNITS 

GPTION NUMBER? (0-6) 

The user note.; that fuel and air use the same units for temperature. Again, 
order Is unimportant, but the user chooses to work top-down through the menu. 
Next In the menu Is fuel temperature selection. 

0 

ENTER TEMPERATURE OF FUEL MIXTURE, DEG- CELSIUS 
The user responds with the desired fuel temperature. 


125.00 


“•CURRENT DATA STATUS*** 


FUELS AND 
FUEL 
C3H8 

FUEL TEMP 
AIR TEMP 
PRESSURE 
A-F RATIO 


FUEL AMOUNTS : 

RELATIVE MOLES AMONG FUELS 
1 .0000000 

= 1 . 2500E+02 DEG- CELSIUS 

= O.OOOOE-Ol DEG- CELSIUS 
= 0.0000E-01 ATMOSPHERES 
= U.OOOOE-Oi FUEL-AIR EQUIVALENCE RATIO 


SELECT AN OPTION: 

0- RUN WITH THIS DATA 

1- CHANGE FUELS OR FUEL AMOUNTS 

2- CHANGE FUEL TEMPERATURE 

3- CHANGE AIR TEMPERATURE 

4 - CHANGE TEMPERATURE UNITS 

5- CHANGE PRESSURE OR PRESSURE UNITS 

6- CHANGE AIR-FUEL RATIO OR RATIO UNITS 


OPTION NUMBER? (0-6) 

m 


ENTER TEMPERATURE OF AIR, DEG- CELSIUS 

» 

Now the desired Input air temperature Is specified. 



“‘CURRENT DATA STATUS*** 

FUELS AND FUEL AMOUNTS : 

FUEL RELATIVE MOLES AMONG FUELS 

C3H8 1.0000000 
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¥ 


FUEL TEMP - 1 .2500E*02 DEG- CELSIUS 

AIR TEMP = 3.7500E+02 DEG- CELSIUS 

PRESSURE = O.OUOOE- 01 ATMOSPHERES 

A-F RATIO = O.OOOOE-Ol FUEL-AIR EQUIVALENCE RATIO 

SELECT AN OPTION: 

0- RUN WITH THIS DATA 

1- CHANGE FUELS OR FUEL AMOUNTS 

2- CHANGE FUEL TEMPERATURE 

3- CHANGE AIR TEMPERATURE 

4- CHANGE TEMPERATURE UNITS 

5- CHANGE PRESSURE OR PRESSURE UNITS 

6- CHANGE AIR-FUEL RATIO OR RATIO UNITS 

OPTION NUMBER? (0-6) 

Still working arbitrarily top-down, only options 5 and 6 remain to be 
exercised before the dataset Is complete. 

a 

SELECT AN OPTION: 

T- SAME UNITS, NEW PRESSURE 

2- CHANGE UNITS TO KILO PASCALS 

3- CHANGE UNITS TO P.S I. 

4- CHANGE UNITS TO A’ .JSPHERES 

5- CHANGE UNITS TO BAR 

OPTION? (1-5) 

Note that entering either a "1" or a "4" would have the same effect here. 

□ 

ENTER PRESSURE. IN ATMOSPHERES 

am 

***CURRENT DATA STATUS*** 

FUELS AND FUEL AMOUNTS : 

FUEL RELATIVE MOLES AMONG FUELS 

C3H8 1.0000000 

FUEL TEMP = 1 .2500E+02 DEG- CELSIUS 

AIR TEMP = 3.7500E+02 DEG- CELSIUS 

PRESSURE = 4 . 5000E+00 ATMOSPHERES 
A-F RA1I0 = O.OOOOE-O! FUEL-AIR EQUIVALENCE RATIO 

SELECT AN OPTION: 

0- RUN WITH THIS DATA 
T- CHANGE FUELS OR FUEL AMOUNTS 

2- CHANGE FUEL TEMPERATURE 

3- CHANGE AIR TEMPERATURE 

4- CHANGE TEMPERATURE UNITS 

5- CHANGE PRESSURE OR PRESSURE UNITS 

6- CHANGE AIR-FUEL RATIO OR RATIO UNITS 
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) 


OPTION NUMBER? (0-6) 


The only null data left Is the a 4 r-fuel ratio. 

0 

SELEC1 AN OPTION: 

1- SAME UN11S, NEW RAUO 

2- CHANGE UNITS TO KG AIR/KG FUEL 

3- CHANGE UNI 1 S TO KG FUEL/KG AIR 

4- CHANGE UNITS TO MOLS AIR/MOLS FUEL 

5- CHANGE UNITS TO MOLS FUEL/MOLS AIR 

6- CHANGE UNITS TO FUEL-AIR EQUIVALENCE RATIO 

OPTION? (1-6) 

Note that either a "1" or a "6" will have the same effect here. 

0 

ENTER FUEL-AIR EQUIVALENCE RATIO 

EH 

‘“CURRENT DATA STATUS“* 

FUELS ANO FUEL AMOUNTS : 

FUEL RELA1IVE MOLES AMONG FUELS 

C3H8 1.0000000 

FUEL TEMP = 1 . 2500E+02 DEG- CELSIUS 

AIR TEMP = 3 . 7500E+02 DEG- CELSIUS 

PRESSURE = 4 . 5000E+00 ATMOSPHERES 

A-F RATIO = 1 .1000E+00 FUEL-AIR EQUIVALENCE RATIO 

SELECT AN OPTION: 

0- RUN W11H THIS DA 1 A 

1- CHANGE FUELS OR FUEL AMOUNTS 

2- CHANGE FUEL TEMPERATURE 

3- CHANGE AIR TEMPERATURE 

4- CHANGE TEMPERATURE UNITS 

5- CHANGE PRESSURE OR PRESSURE UNITS 

6- CHANGE AIR- FUEL RATIO OR RA1 10 UNI IS 

OPTION NUMBER? (0-6) 

The user has now specified all Input conditions of the problem. 

0 

RUN WITH THIS DATA: ARE YOU SURE? (Y/N) 

The user has double- checked all the Items on the menu, so as to avoid wasting 
computer time on bad Input data. 

0 


ij 
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EQLBRM executes the dataset and presents the solution: 


i 


OUTPUT 

****** 


*** ADIABATIC EQUILIBRIUM AT C0NS1ANT PRESSURE AND ENTHALPY *** 

PRESSURE= 4.500E+00 ATMOSPHERES 
TEMPERATURE = 2007.09 DEGREES CELSIUS 

FUEL-AIR EQUIVALENCE RATIO = T.TOO 


SPECIES 

MOLE 

MOLE 

MASS 

EMISSION 

NAME 

NUMBER 

FRACTION 

FRACTION 

INDEX 

("IM 

(MOL-I/G-X) 

(MOL-I/MOL-X) 

(G-I/G-X) 

(G-I/KG-FUEL) 

CO 

1 .030E-03 

2.993E-02 

3 . 026E-02 

3 .026E+01 

C02 

3. 109E-03 

9 . 444E-02 

1 .500E-01 

1 .500E+02 

H 

1 . 407E-05 

3 . 899E-C4 

1 .418E-05 

1 .418E-02 

H2 

3 . 274E-04 

9.071E-03 

6 . 600E-04 

6.600E-01 

H20 

5.632E-03 

1 .560E-01 

1 .015E-01 

1 .01 5E*02 

N 

4.476E-10 

1.240E-08 

6.269E-09 

6.269E-06 

NO 

1 .784E-05 

4 . 944E-04 

5.354E-04 

5 . 354E-01 

N2 

2.556E-02 

7 . 083E-01 

7.161E-01 

7 . 1 61 E*02 

0 

1 .125E-06 

3.11 8E-05 

1 .800E-05 

1 .800E-02 

OH 

3.805E-05 

1 .054E-03 

6.472E-04 

6 . 472E-01 

02 

8.201E-06 

2 . 272E-04 

2 . 624 E - 04 

2.624E-01 


DO YOU WANT A HARD COPY OF YOUR INPUT AND OUTPUT? (Y/N) 

The user elects to make a hardcopy of both the Input and output datasets. 

0 

DATA AND RESULTS ARE BEING APPENDED TO 
FILE "EQLOUT.LST" IN YOUR DIRECTORY. 

The hard ropy output Is directed to logical unit 2. System- specif 1c commands 
must be executed at run time (or added to the source code) to assign this unit 
to the file “EQLOUT.LST," or directly to an output device. 

2.3 Debugging 

If Improper operation of EQLBRM Is suspected, a debugging feature can be 
brought up by "un-commentlng" the appropriate source code statements near the 
top of the main routine. With this feature restored, an Interactive choice of 
debugging level (1-5) can be selected, which can be used (together with the 
source code listing) to discover Implementation or source code transcription 
errors . 

2.4 FOR' RAN Language and Availability 

Program EQLBRM was developed on a Digital Equipment Corporation POP 11/23 
computer. The source code Is In 1966 ANSI-standard FORTRAN IV, and runs 
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without modification (except for system specific 1/0 statements In the main 
routine KEYBRO) on VAX and IBM computers. 

Source codes and the required data file EQLBRM.DA1 are available from 
COSMIC, University of Georgia, Athens, Georgia 30602. 

3. THERMOCHEMICAL D.aI A 

lhe thermochemical data for elements and for nonelemental molecular species 
are taken directly from the NASA element and thermochemical data files (ref. 1). 

In the EQLBRM.DA1 data file, the ELEMENTS code word and elemental data 
must precede the THERMO code word and thermochemical data. The format for 
input data are summarized in tables I and II: 
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TABLE I. - ELEMENTS INPUT RECORDS 


Order 


Contents 

Format 

Record 

Columns 

First 

ELEMENTS 

3A4 

1 to 8 

Any 

One record each distinct 
element present In the chemical 
system. Each record contains: 




1) 

Atomic symbol of element-- 
must agree with that used 
In THERMO data. 

A2 

1 to 2 


2) 

Atomic weight of the element 

F10.6 

10 to 19 


3) 

Valence or oxidation state 
of the element (positive, 
negative or zero) . 

F10.6 

20 to 29 



See Note A 



Last 


Blank record card 


all 


Note A: EQLBRH is limited to ten distinct elements. 



TABLE II. - THERMO INPUT RECORDS 


Order 


Contents 

Format 

Card 

Columns 

First 

THERMO 

3A4 

1 to 6 

See 

Sets 

of four records In sequence for each 



Note B 

species In the chemical system. The 




record formats for each set are, In order: 



; First 

1) 

Molecular symbol or name of species 

3A4 

1 to 12 

| In set 

2) 

Date 

2A3 

19 to 24 


3) 

Atomic symbols and formula 

4(A2,F3.0) 

1 


4) 

Phase (gas only, letter G) 

A1 

45 

I 

5) 

Temperature 

2F10.3 

46 to 65 


6) 

Integer 1 

115 

80 

Second 

1) 

Coefficients ( , 1 = 1,5) for upper 



In set 


temperature range. See Note A. 

5E15.8 

1 to 75 


2) 

Integer 2 

15 

80 

1 Third 

1) 

Coefficients Z 5 and Z 7 for upper 



1 In set 


temperature range, and z-| , Z 2 and 





Z 3 for lower. See Note A. 

5E15.8 

1 to 75 


2) 

Integer 3 

15 

80 

! Fourth 

1) 

Coefficients (z^, 1 = 4,7) for low 



In set 


temperature Interval. See Note A. 

4E15.8 

1 to 60 

| 

2) 

Integer 4 

120 

80 

1 


See Note C 

120 

80 

Last 


Blank record 




Note A: Coefficients (zj , i = 1,7) correspond to the following equations for constant- 

pressure specific heat capacity, enthalpy, and entropy: 


C 

_£ 

R 


2 3 

+ z T + z.T * z T 
2 3 4 


V 



ORIGINAL PAGE 19 
OF POOR QUALITY 


s° 

R = Z 1 


log T 


V 


Z 3 T 


Z 4 T 


Z 5 T 

+ T~ * Z 7 


The individual species enthalpy h include the standard-state enthalpy of formation and 
the sensible enthalpy. The individual species entropy s° is for a pressure of one 
atmosphere. 

Note 8 : The first two sets must be molecular nitrogen N 2 and molecular oxygen 0 2 , the 

only constituents of air assumed in EOLBRM. All polyatomic fuel data-sets must reside at 
the end of the THERMO data file. 

Note C: EQIBRM is limited to 80 distinct atomic and/or molecular species, including those 

entered interactively at run time. 
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